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Abstract 

We study the most general case of spherically symmetric vacuum 
solutions in the framework of the Covariant Horava Lifshitz Gravity, 
for an action that includes all possible higher order terms in curva- 
ture which are compatible with power-counting normalizability re- 
quirement. We find that solutions can be separated into two main 
classes: (i) solutions with nonzero radial shift function, and (ii) solu- 
tions with zero radial shift function. In the case (ii), spherically sym- 
metric solutions are consistent with observations if we adopt the view 
of Horava and Melby-Tomson [II], according to which the auxiliary 
field A can be considered as a part of an effective general relativistic 
metric, which is valid only in the IR limit. On the other hand, in the 
case (i), consistency with observations implies that the field A should 
be independent of the spacetime geometry, as the Newtonian poten- 
tial arises from the nonzero radial shift function. Also, our aim in this 
paper is to discuss and compare these two alternative but different 
assumptions for the auxiliary field A. 
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1 Introduction 



A recent higher order space derivative model for Gravity was formulated 
by Horava [lj. This model is power-counting renormalizable and serves as 
an ultraviolet (UV) completion of General Relativity (GR). This scenario 
is based on an anisotropy between space and time coordinates, which is 
expressed via the scalings t — > b z t and x — > bx, where z is a dynamical 
critical exponent. It is worth noting that Horava-Lifshitz (HL) Gravity has 
stimulated an extended research on Cosmology and black hole solutions, see 
for example [2]- [27], and we note that quantum field theory models in flat 
anisotropic space-time were also developed, see for example [2B]-[23] and 
references therein. 

In HL Gravity, the four- dimensional diffeomorphism invariance of GR is 
sacrificed in order to achieve power-counting renormalizability. Although HL 
Gravity violates local Lorentz invariance in the UV, GR is expected to be 
recovered in the infrared (IR) limit. This implies a very special renormaliza- 
tion group flow for the couplings of the model, in particular it is expected 
that the coupling A in the extrinsic curvature term of the action has the 
behavior A — > 1, i.e. that it flows towards its GR value. But there is no 
theoretical study supporting this specific behavior. In addition, there are 
several other potential inconsistencies in HL Gravity which have been dis- 
cussed (see for example [M]-[in] and references therein). More specifically, 
the breaking of 4D diffeomorphism invariance introduces an additional scalar 
mode which may lead to strong coupling problems or instabilities, and in this 
way prevents HL Gravity from fully reproducing GR in the IR limit. 

In Ref. jH] a new Covariant HL Gravity is formulated by Horava and 
Melby-Thompson (HM), which includes two additional nondynamical fields A 
and u, together with a new U(l) symmetry. In this model the U(l) symmetry 
eliminates the extra scalar mode curing the strong coupling problems in the 
IR limit. Note, that in what follows we adopt the HM action of Ref. [4T| , 
in which the parameter A is set equal to one (A = 1). However, the U(l) 
symmetry can not force the value of the parameter A to be equal to 1, as an 
action, with the extended gauge symmetry and A ^ 1, has been formulated 
in [12] . Also, HM theory reproduces many features of GR at long distances 
as it is shown initially in Ref.|41j. 

HL Gravity can be separated into two versions which are known as pro- 
jectable and non-projectable. In the projectable version the lapse function 
iV (see Eq. ([2]) below) depends only on the time coordinate, while in the 
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non-projectable version N is a function of both space and time coordinates. 
The Covariant HL Gravity considers the projectable case. 

In this paper we study the most general case of spherically symmetric 
vacuum solutions in the framework of Covariant HL Gravity, for an action 
which includes all possible terms allowed by renormalizability requirement. 
We find that spherically symmetric solutions can be separated into two main 
classes: (i) solutions with nonzero radial shift function, and (ii) solutions 
with zero radial shift function. 

We would like to note that Covariant HL Gravity, as it is formulated 
by Horava and Melby-Tomson (HM) in [H] , incorporates an additional as- 
sumption for the field A. In particular, the field A is assumed as a part of 
an effective general relativistic metric in the IR limit, via the replacement 
N — > N — A/c 2 . Spherically symmetric solutions, in the case (ii), are consis- 
tent with observations if we adopt HM approach for A. On the other hand, 
in the case (i), we see that A should be independent of the lapse function 
in order to achieve consistency with observation. We would like to warn the 
reader, that these two alternative but different views for the field A are exam- 
ined both in this paper, in the framework of spherically symmetric vacuum 
solutions of cases (i) and (ii). 

The paper is organized as follows. In Secj2] we summarize features of the 
Covariant HL Gravity. We consider then in Secj3] the most general ansatz 
for spherically symmetric solutions, including a nonzero radial shift function, 
and we derive the equations of motion and the corresponding constraints. In 
SecJUwe present the solutions for the two situations (%) and (ii), and finally 
Secj5] contains our conclusions. In Appendix A are presented some details on 
the geometrical interpretation of the auxiliary field A. In Appendix B, we 
present the minimal substitution approach of Ref. [12] which describes how 
mater couple with the auxiliary fields A and v. 

1 For the construction of the HL action, the so called "detailed balance principle" has 
been proposed jT]. The main advantage of this approach is the restriction of the large 
number of arbitrary couplings that appear in the action of the model. However, a more 
general way for constructing the action would be to include all possible operators which are 
compatible with the renormalizability [21113]; this implies that all operators with dimension 
less or equal to six are allowed. 
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2 Covariant Horava-Lifshitz Gravity 



In this section we introduce the notation for the Covariant HL Gravity in the 
case of three spatial dimensions (d = 3), and we discuss the additional U(l) 
symmetry of the model, as well as the role of the auxiliary nondynamical 
fields A and v which are included in the action. 

2.1 The action 

This model, as the original HL Gravity, is characterized by an anisotropy 
between space and time dimensions 

[t] = -z, [x] = -1 , (1) 

where z is an integer dynamical exponent. The action of the model is struc- 
tured by a set of five fields: N(t), Ni(x,t), gij(x,t), A(x,t) and u(x,t) 
(i = 1,2,3). Note that N(t), Ni(x,t), gij(x,t) are the standard fields that 
appear in the Arnowitt, Deser and Misner (ADM) form of the space-time 
metric 

ds 2 = -c 2 N 2 dt 2 + gij (dx i + Wdt) (dx j + N j dt) , (2) 

where cis the velocity of light, with dimension [c] = z — 1, and spatial compo- 
nents dx l /dt (i = 1,2,3). In addition, iV and iVj are the "lapse" and "shift" 
functions which are used in general relativity in order to split the space-time 
dimensions, and is the spatial metric of signature (+,+,+). Note that 
here we are interested for the projectable version of the model which implies 
that the lapse function N(t) depends only on the time parameter. For the 
dimensions of "lapse" and "shift" functions we obtain 

[N] =0, [Ni] = z-1. (3) 

The auxiliary fields A(x, t) (potential) and v(x, t) (prepotentiald ) are nondy- 
namical fields which have to satisfy constraint equations. As we will see sub- 
sequently the existence of these fields is necessary in order to achieve invari- 
ance of the action under the extended Gauge symmetry: U (1) x Diff(Ai,J r ). 
In addition, the dimensions of these fields are 

[A] =2z-2, [u] = z-2. (4) 

2 The terminology "potential" for A(x,t) and " prepotential" for u(x, t) has been intro- 
duced in the original work of |41j 
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The full action of the model is formulated as 

S = A J dtd 3 Xy/g^N [KijK ij -K 2 -V + vQ ij {2Kij + V^Vj 

-A(fl-2fi)}, (5) 

in which d is the spatial dimension (D = d+ 1 = 4), k 2 is an overall coupling 
constant with dimension [k 2 ] = z — d, and the extrinsic curvature is 

Kij = ^ - V^- - VjiVJ , z, j = 1, 2, 3 (6) 

where the symbol G 1 - 7 is defined as 

Q ij = R ij - X -Rg 13 + Qg ij . (7) 

Note that this choice for z = 3 is an immediate consequence of power count- 
ing renormalizability request. In particular, the coupling k 2 in the above 
action has dimension [k 2 ] = z — 3, hence, if z — 3 the HL Gravity model 
is renormalizable, for z > 3 is super-renormalizable and for < z < 3 is 
non-renormalizable . 

For the construction of the potential term V we will not follow the stan- 
dard detailed balance principle, but we will use the more general approach 
[21 H3], according to which the potential term is constructed by including 
all possible renormalizable operators (relevant and marginal) □, that have 
dimension smaller than or equal to six, hence we write 

V = V IR + V R 2 + V R 3 + V AR 2 (8) 

where 

V IR = -c 2 (R -2 A) (9) 

V R 2 = - ai R 2 - a 2 R ij Rij 

V R 3 = — [3iR 3 — (3 2 RR lj Rij — PsRlRjRl 

V AR 2 = -f3 4 RV 2 R - faViRjkVRi" 

The dimensions of the various terms in the Lagrangian are 

[R] = 2, [R 2 ] = 4, [R 3 ] = [AR 2 ] = 6, (10) 



3 We have ignored terms which violate parity, see also [43] . 
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where the symbol A is defined as A = gj J 'V l V- J (i = 1, 2, 3). In addition, we 
have used the notation R, Rij and Rijki for the Ricci scalar, the Ricci and the 
Riemann tensors = 1,2,3), which correspond to the spatial 3D metric 
gij. Note that the term R t i H R i jki does not appear in Vr?, as the Weyl tensor 
in three dimensions automatically vanishes. In addition, c 2 , the couplings q.{ 
{% = 1,2), and (3j (j = 1, • ■ ■ , 5), have dimensions 

[ c 2 ]=4 , h = 2 , [/y = 0. (ii) 

Finally, we would like to note that the potential terms, of Eqs. (jSJ) and (jUJ) 
above, has been considered previously for covariant HL gravity, in the case 
of cosmology, by the authors of Ref . [H] . 



2.2 The extended U(l) X Diff(M,J r ) Gauge symmetry 

The main motivation for considering models with an anisotropy between 
space and time dimensions type is the construction of a power-counting renor- 
malizable Gravity model. However, in order to achieve normalizability, and 
simultaneously keep the time derivatives up to second order, we have to sacri- 
fice the standard 4D diffeomorphism invariance of General Relativity, which 
is now restricted to 

6t = f(t), 5x l = at,^) (12) 

which is a foliation preserving diffeomorphism, Di f f (Ai , J 7 ) , where M. is 
the spacetime manifold, provided with a preferred foliation structure J 7 . In 
particular the fields N, iVj, transforms as: 

Sgij = d£ h g ik + d^ k g ik + £ k d k g i3 + fg i3 (13) 
SNt = dgNj + gdjNi + ggij + fNi + fNi 
SN = gdjN + fN + fN 

However, the action of Eq. (jSJ) has an additional symmetry, in particular 
it remains invariant under a U(l) Gauge symmetry, according to which the 
fields of the model transform as 

5 a N = 

S a 9ij = 

8 a Ni(x,t) = NVia 
5 a A(x,t) = a-N { Via 

5 a v = a (14) 
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where a is an arbitrary spacetime function. Accordingly, the full symmetry of 
the action of Eq. ^ is the extended Gauge symmetry: U{\) x Diff(M,J c ). 

The extrinsic curvature term in the action of Eq. ([5]), if we assume only 
Dif f{Ai,J z ) symmetry, could include an additional running coupling con- 
stant A appearing as: 

C K = K^K t] - \K\ (15) 

In order to achieve agreement with General Relativity, we expect that RG 
flow in the IR leads the coupling A to unity, but the exact mechanism for 
this remains unknown. There was a hope, in Ref. [H], that the Covariant 
model requires A = 1 due to the U(l) symmetry, hence the above mentioned 
problem, for the flow of A in the IR, does not exist. On the other hand, in Ref. 
[12] , it was shown that the action of Eq. fl5]) can be written in an invariant 
form under U(l) symmetry, for an arbitrary value of A, which implies that 
U(l) symmetry can not fix the coupling A (see Appendix B in the present 
paper). 

Note that according to HM approach an infinitesimal U(l) transformation 
on the fields A, u, N, N l , in the nonrelativistic limit, is equivalent to an in- 
finitesimal diffeomorphism involving the time coordinate. As a consequence, 
the symmetry U(l)x Diff(A^, J 7 , ) can be seen approximately as a Diff(M,4) 
symmetry of Standard General Relativity, see Ref. [H] and the discussion in 
our Appendix A. 

2.3 Infrared limit 

Although the additional U(l) symmetry eliminates the extra degree of free- 
dom, the IR limit of the HM theory, obtained after neglecting higher order in 
spatial curvature terms, does not coincide with General Relativity. Indeed, 
the action is then 

S = [ dx°d 3 x^j{N [K l3 K %j -K 2 + R + v&\2K i3 + V.V^)] 

lo7rGr J L 

-A(R-2Q)\ (16) 

where the time-like coordinate Xq is defined as Xq = ct, and the fields (in 
the above action) are rescaled according to Sec. 5 in Ref. |41j . However, as 
discussed in [41] , Covariant HL Gravity reproduces many features of general 
relativity for long distances. 
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3 Spherically symmetric solutions in Covari- 
ant HL Gravity 

3.1 The metric 

The starting point is the action describing Gravity with anisotropic 
scaling: 

S = A J dtd 3 Xy/g^N [KijK ij - K 2 - V + vQ ij {2K ij + ViV»] 

-A(fl-2fl)}, (17) 

We consider the most general static spherically symmetric metric, of the 
form: 

ds 2 = -c 2 N 2 dt 2 + — — (dr + n(r)dt) 2 + r 2 (sin 2 QdQ 2 + d<f) 2 ), (18) 

fin 

where n(r) = N r (r) is the radial component of shift functions, and iV r = 
n(r)/f(r) since g„ = l//(r). 

3.2 Constraints 

The auxiliary fields v and A lead to the following constraints 

• The variation of S with respect to A gives R — 2Q = 0, or equivalently 

J2 = -^(r/' + /-l) = 2n , (19) 
which imposes that the function f(r) in the metric ffl8|) is 

. . fir, IB 

f(r) = 1 - -r 2 , (20) 

where 5 is a constant of integration. As we will see, in the following 
sections, this constant of integration B is not interpreted necessary as 
the mass of the spherically compact object. Moreover, in what follows, 
we will set Q = 0, as we are looking for asymptotically flat solutions. 
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• The variation of S with respect to v gives 

(-)''V,V ; // • (-)'•'' 0. (21) 

In what follows we will assume the Gauge fixing v = 0, then the above 
constraint gives Q lJ Kij = 0, which is satisfied for spherically symmetric 
solutions as we will see in the following sections. 



3.3 Equations of motion 

We therefore start with the action 



S~ J dtd 3 x^{N (T-V)- AR], T = K ij K ij -K 2 1 (22) 

in which we have considered the Gauge fixing v = 0, and we have set Q = 0. 
The Lagrangian which corresponds to the action (1221) . has the form 



2 



where 



£ = C T + C V -^==AR, (23) 

my , 

£t = r — T{r,n,n,f,f), 

Vf( r ) 

N(t)r t n w 

Cv = Trr^VviJiJ »/ >/ )> 

v/( r ) 



T = -j^r 2 (n 2 f + 2rnfn' + rn 2 ff) 
r 2 . „ 2 



r-^AR = - 7 ==A(r)(rf' + f-l). (24) 

The potential term of the Lagrangian, for zero cosmological term (A = 0) 
(as we are interested for solutions which are asymptotically flat), is taken as 

V= - c 2 R- ai R 2 -a 2 R ij Rij - (3iR 3 - f3 2 RR ij Rij (25) 
- p 3 K j R? k R k i - faRV 2 R - p B ViR jk V i R ik , 
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in which 



R = --Arf-l + f) 



R ij ^j = ^i[3(r/') 2 + 4rr(/-l)+4(/-l) 2 ] 

R 3 R\R\ = -^[5(r/') 3 + 6(r/') 2 (/-l) + 12r/'(/-l) 2 + 8(/-l) 3 ] 

RV 2 R = 2(r// ~ 1 + /) x 

[2r 3 /"7 + r 3 /7' + 2r 2 /7 + 2(1 - 3/)r/' - 4/(1 - /)] 
V^ViP* = l E [3(r 2 f) 2 -2r 3 ff" + Sr 2 f"(l-f) + 

+5(r/') 2 + 16(1 -f)rf' + 24(1 - /) 2 ]. 
The Euler equation for n gives 

d (dC T \ dC T 



dr V <9n' / dn 
from which we obtain 

f(r)n(r) = 0, (27) 

such that necessarily either n(r) = or /(r) is constant. The Euler equation 
for / is 

n=0 v J 7 

where /( n ) = d n f/dr n , hence we obtain 

*+sH) + ^-°* (29) 

where a prime denotes a derivative with respect to r and the differential 
operator O is 

4/ 2r ^ ( 1} dr» W70/W,/- (30) 



9 



After same algebra, for f(r) — 1 — and N — 1, we obtain 

£>5 = - - ' _ (31) 
r(r - 25) y ' 

OR 2 = 

B 2 



ORR^Rij 



OR i j R j k R k 



2r 4 (r - 25) 
OR 3 = 

65 2 (115 - 6r) 

r 7 (r -25) 
35 2 (505 -27r) 



2r 7 (r - 25) 

(95V 2 5 = 

fc 352(405 - 21r) 



OV t R jk V l W 



2r 7 (r - 25) 



Finally, the variation of the action with respect to N(t) gives the so called 
Hamiltonian constraint 



r 2 



dr^==(T + V) = 0, (32) 

and, using a time redefinition, the lapse function N(t) is set equal to unity 
(N(t) = 1). 



4 Analytic solutions 

We have to satisfy the following two equations of motion 

/' (r)n(r) = (33) 

2r\fJ 



(34) 



and the Hamiltonian constraint 

roo 2 

/ dr-=(T + V) = 0, (35) 
Jo \/.f(r) 
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while the constraint of the spatial curvature (R = 0) gives 

f(r) = 1 - — . (36) 
r 

The constraint for v (Eq. (j2ip ) is verified for ^ = 0, what we will assume in 
what follows. 

4.1 Nonzero radial shift function: n(r) ^ 

If the shift function does not vanish, we see from the constraint (133]) that / 
must be constant, and thus B = 0, such that 

/ = 1. (37) 

As a consequence, V = (see Eqs. (13~T]) ). and Eq. (}3"4"]) can be written as 

rA' + 2{rn 2 )' = 0, (38) 

while Eq. gives the Hamiltonian constraint: 

dr(rn 2 Y = . (39) 

From Eqs. fl38|) and f )39|) . we obtain then 

i 

dr rA'(r) = . (40) 



o 





The above integral is convergent only if A(r) has the following large and 
short distance asymptotic behavior 

• A(r) ~ Coo + Ca r~ b for r -)■ +oo, with 6 > 1; 

• A(r) ~ C + C A r~ a for r ->• 0, with a < 1; 

in which Coo, C^, Co and C^ are arbitrary constants. An integration by 
parts, of Eq. ( l4"Uj) . leads to 

oo poo 

drrA'ir) = / drr (A(r) — Coo)' 
n Jo 

oo 



oo 



[r(A(r)-Coo)]o - / ^(A(r)-Coo) 

rfr (A(r) - Coo) (41) 
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Note, that the fields A and A + const are equivalent in the sense that they 
give the same n 2 from Eq. ( 13 8 p . hence without loss of generality we can set 



Coo = hm A(r) = (42) 

r— >oo 

In this case the Hamiltonian constraint of Eq. (j4*0|) is equivalent to 

drA(r) = (43) 



In what follows we consider separately two cases 1) A(r) = which is the 
minimal choice, 2) A(r) is a function which satisfy the Hamiltonian constraint 
of Eq. (j43|) (or the equivalent equation (j40|) ). 

4.1.1 First case: A=0 

When A = 0, the constraint of Eq. ([43]) is satisfied, hence from Eq. (1381) we 
obtain 

n(r) = ±y^, C M = 2GMc 2 (44) 

where Cm is a constant of integration. Thus the metric of Eq. ( 1T81) can be 
written as 



rfs 2 = _ c 2 dt 2 + \ dr± f^M^ dt \ + r 2( sin 2 0d0 2 + ^ 2)> (45) 



which is the Schwarzschild solution in Painleve-Gullstrand coordinates, see 
for example Ref. [1] and references therein. Note that in this case, the 
Newtonian potential <p(r) (goo = 1 + 20(r)) is proportional to the square of 
the radial shift function, according to the equation 



n\r) _GM 



2c 2 r 

However, we observe that the above expression for the Newtonian potential 
is not U(l) invariant. In particular it is the Newton Law for particular 
Gauge choice v = 0. We can correct this situation if we take into account 
the coupling between shift functions iVj and matter. In Ref. [42] (see also 
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Appendix B) we see that iVj couple with matter in the U(l) invariant form 
Ni — NViV, such that the Newtonian potential should be modified as 



<f>(r) 



(n(r) — V r z/(r)) 2 
2^ 



CM 



r 



(47) 



The above expression is U(l) invariant, and for the Gauge choice v = we 
recover Eq. ( 146]) . 

Finally, in the case of the initial HL Gravity model with projectability 
condition, spherically symmetric solutions have been studied in Ref. [4j. The 
solution we present in this section, for Covariant HL Gravity, is also a solution 
of HL Gravity without the U(l) symmetry, but these two models have not 
the same full spectrum of solutions. Also, spherically symmetric solutions 
with nonzero energy momentum tensor (stars), in the case of projectability 
condition, have been studied in Ref. [15] , and spherically symmetric star 
solutions are studied in Ref. [17]: the main conclusion is that a spherically- 
symmetric star should include a time-dependent region near the center. 

4.1.2 Second case: f °° drA{r) = 

In this situation, the most general solution of Eq. fl3"5|) is 



where Cm is a constant, and we will consider separately the two cases: (a) 
C M = and (b) C M + 0. 

(a) We have for C M = 



Since A{r) behaves as A{r) ~ — Cat b for large distances, Eq. (149]) gives a 
modified Newtonian potential when r — > oo 





(48) 




(49) 



4>{r) 





(50) 
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in which ^ 

C * = 4c 2 (6-l)' (51) 

where Ca > is dimensionful, in order for the potential to have the correct 
dimensionality, and is interpreted as the mass of the compact object. 
Also, because n 2 > 0, we have the additional constraint: 

IT 

\ A(p)dp>rA(r), (52) 
Jo 

and to show a function A(r) satisfying the above constraints can be found, 
we give here two examples. 

A first example of function A(r), which satisfies the condition ( 14T1) . is 

Mr) = -^(l-$ i ), b>l (53) 

where 



sin (? - Sf) 
sm If J 

The condition fl52|) is then satisfied for all r only if Ca > 0, such that the 
potential fl5U|) is negative, as expected. Note that this function is singular for 
r = 0, and another example of function A(r) which is regular at the origin 
and which satisfies the constraints, is: 

Mr) = - T ff Eri (r-l2), b>l (55) 

where 

sin (izt) 

72 = — (56) 

Although we have not check it, by choosing b appropriately closely to unity 
b ~ 1 it may be possible the above two solutions to pass solar system tests. 
However, the purpose of this paper is not to set constraints on the parameter 
b, hence the topic of constructing solutions which can satisfy solar system 
tests is left for future investigations. Also, note that the leading order of 
large distance behavior for the radial shift function n(r) (or A(r)) is fixed 
by the requirement of recovery of Newton Law in the large distance limit, 
but for small distances there is an ambiguity in the exact shape of n(r) (or 
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A(r)), as we see from the above two examples. 



(b) For Cm ^ we have 

n 2 (r) = ^--A(r)--^ dp A(p), (57) 

which corresponds to a qualitatively different situation from the one where 
Cm = 0. Here the constant Cm is proportional to the mass of the spherical 
compact object Cm = Cm = 2GMc 2 , and the auxiliary field A determines 
the subleading behavior in the asymptotic expansion of n 2 for large r. As in 
the previous case (Cm = 0), we can choose suitably the field A(r) in order 
to satisfy the Hamiltonian constraint of Eq. (HIT) , and the restriction n 2 > 0. 
It is explained in [IB] that solar system tests requires a large distance asymp- 
totic behavior for A of the form: 

C A 

A(r) ~ —j- + • ■ ■ , when r — > oo, (58) 

where the exponent satisfies b > 3, and the dots represent higher order pow- 
ers of 1/r. We find here that the following choice for A(r) has the required 
asymptotic behavior, and satisfies the corresponding constraints (the Hamil- 
tonian constraint of Eq. (14ip and the requirement n 2 > 0), 

Mr) = ^(l-Tar 62 ) , (59) 
with 3 < bi, 3 < bi — b 2 , and — 1 < b 2 ^ 

where is a dimensionful constant and 

sin ( + 2k) 
73 = ^7 y . (60) 



sin | Vi 



Hence solar system tests do not necessarily impose A = 0, and the choice 
for the function A remains an open question. Furthermore, according to the 
above examples we see that there is a freedom in the choice of A, which cor- 
respond to different spherically symmetric solutions, in contrast to standard 
General Relativity for which spherical symmetry (in the absent of matter) 
leads to Schwarzschild geometry. 
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4.2 Zero radial shift function: n(r) = 

In this situation, /(r) = 1 — 2£>/r, with B ^ 0, and the evolution equation 
fl29|) for A gives 



eL4 c 2 5 a 2 5 2 



dr r(r - 2B) r(r - 25) 2r 4 (r - 2B) 

6/3 2 5 2 /6 UB\ 



(61) 



(r-2B) \r 



mB 2 



20-25 



.6 




27 


505 






21 


405 


j -6 





The solution of this equation is 

A{r) = c 2 + A VT^-^- 2 {-2 + 2x + x 2 + x 3 ) + 6 -^x e 
fa ( 4 , 4 - , 2 _ 2 , 2 3 , 5 4 , x 5 



'■ >> h (A + \ x + 2 2 + 2 3 + 1^4 + ± + 20x 6 j ^ ( g 2) 



US 4 V 7 7 7 7 14 2 

where x = B/r and v4 is a constant of integration. If we set B = GM Q, the 
solution of Eq. (162]) can be written as an expansion in x = M/r (from now 
on it is convenient to set G = 1) 

m \ 2 . A «2 4/3 3 12/3 5 

AW = C + ^ + 5M* " 77M* + 77M* (63) 

( x 2 x 3 \ ( . a 2 4/3 3 12/3 5 ^ 
+ x + — + — -A - — — + 



2 2 J \ u 5M 2 77M 4 77M 4 J 

4 f 5^ 15& \ 

! 154M 4 154M 4 J 



+x 5 



7Ao ^3_ _ 3/3 5 \ 
8 22M 4 22M 4 J 



4 In the case of solution with zero shift function the constant of integration B is pro- 
portional to the mass of the compact object. 
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1 8 M 4 11M 4 11M 4 y 
+0(x 7 ). 

in which the constant of integration A is chosen to vanish the constant term 

C +Ao+ 5A^-lli^ + 7W =0 ' (64) 
Now A(r) can be split to an IR and an UV part according to the equation: 

A(r) = A IR (r) + A uv (r) (65) 

where we have set 



A IR (r) = c 2 (l-VT^2x~) (66) 
A f \ ( " 2 4 & 12 & \ r, 7T 

"2 i n , n __ , _ 2 , _3\ , 6/^2 6 



(-2 + 2x + x 2 + x 3 ) + 



10M 2 v ' M 4 

ft / 4 4 2 2 2 3 5 4 x 5 nr e 
H — ^-T77 K -x + -x 2 + -x 3 H x 4 H h 75x e 



11M 4 V 7 7 7 7 14 2 

3/3 5 / 4 4 2 2 2 „ 5 , x 5 , . , _ 

-- + =x + -x 2 + -x 3 + — x 4 + — + 20x 6 (67) 



11M 4 V 7 7 7 7 11 2 
4.2.1 The potential interpretation of A 

According to the original formulation of covariant HL gravity by Horava 
Melby (HM) in |4T] the field A has a particular role in the IR limit: more 
specifically it is promoted as a part of an effective general relativistic metric 
via the replacement 

N -> N - ^ . (68) 
c 2 

Hence the spacetime geometry is determined effectively by the following met- 



ric 



ds 2 = -c 2 (^N 2 - N N i+J A iR N ^j dt 2 + 2N i dx i dt + gij dx l dx j , (69) 
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which is realized in the nonrelativistic limit c — > oo, by dropping higher order 
terms in 1/c 2 , as it is shown initially in [JT] (for details see also Appendix 
A). 

In the above effective metric of Eq. fl69|) we have included only the Ajr 
part of A which is consistent with HM approach. On other hand the Auy 
part of A can not be included in this metric, as it is subdominant in the 1/c 2 
expansion. Note that an investigation of the physical role of Ajjy (in HM 
approach) is beyond the scope in this paper. 

It is worth noting that, in the case of solutions with zero shift function 
(Ni = and N = 1), the HM approach allows a potential interpretation for 
the field A. If we take into account that in the IR limit (r — > oo) 

£00 = 1 + 20 + ■ ■ • 

and compare with Eq. (|69|) above, we obtain that the Newtonian potential 
4>{r) is related with A{r) according to the equation 

m = -6up. = -?E+o(m\ (to) 

& r \ r J 

Also the effective metric of Eq. ( |69|) in the IR limit, in the case of solutions 
with zero shift function, can be written as 

ds 2 eff = -c 2 (l - ^\ dt 2 1 dr 2 + rW . (71) 

It is clear that the above metric mimics Schwarzschild geometry in the IR, 
such that the metric of Eq. ( I7ip passes solar system tests, and no restriction 
on Horava couplings is necessary (for details see also [2]). At this point, it 
should be emphasized that the geometric approach for A is an independent 
assumption of HM theory, however in this paper we use it mainly to make 
our solutions (with zero shift function) physically relevant. Finally, the full 
expression of A(r) for r — > oo, is given by the expression 

A(r) M M' 2 M' A 
c 2 r 2r 2 2r 3 

'5 a 2 \ / M 4 7M 5 

1 8 + 8M 2 c 2 ) V~ + ^ 
21 21a 2 6/3 2 147/5 3 lll/3 5 \ M 6 



40M 2 c 2 M 4 c 2 22M 4 c 2 22M 4 c 2 J r 6 
+0(M 7 /r 7 ) . 
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We note that the first two corrections, in M 2 /r 2 and M 3 /r 3 , are independent 
of the Horava-Lifshitz couplings. Also, one can see that corrections of high 
order derivatives vanish in the nonrelativistic limit (c — > oo). 



4.2.2 The Hamiltonian constraint 

In the situation where n — 0, the Hamiltonian constraint reads 

2 

dr ?—V = 0, (73) 
V / 

For the case we consider, f(r) = 1 — 2M/r, the potential V is 

V = ^-^-^--), (74, 

In order to keep the integrand function in Eq. (1731) real, we have to intro- 
duce a spherical gap in space, centered on the black hole and including the 
horizon, assuming that the radial coordinate r has the minimum value L, 
where L > 2M. An equivalent alternative is for example to introduce the 
new coordinate u > 0, defined by r = y/u 2 + L 2 , and to express the whole 
problem in terms of u instead of r. Note that space in this situation is still 
simply connected, such that this non-vanishing length L does not introduce 
topological defects. In this case we have 

dr —=V = QM / 7~7o vU5 + ~T^V ~ 2M ) 

= M~ 1 a 2 C 2 — M~ 3 /3 3 C 3 — M~ 3 /3 5 C 5 , (75) 
where the dimensionless functions C 2 ,C 3 , C 5 are given by the equations 

C( V ) = 4 -- 2 -^ {3 + 2y + 2 ^ (76) 

= |[-|lP (63 + 35! ' + 20! ' 2 

+ I2y 3 + Si, 4 + Si, 5 ) 

C 5 (y) = ^ J 2 (315 + 140^ + 60^ 

+ 24y 3 + 8 i7 4 ), 



3 The authors would like to thank Alex Kehagias for this suggestion. 
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where we have set y = L/M. The Hamiltonian constraint leads to the 
following algebraic equation 

C 2 (y)-kCs(y)-kC 5 (y) = 0, (77) 

in which we have set (3 3 = (3 3 a.2 1 M~ 2 and (3$ = Ih^a^ 1 M~ 2 . We can 
determine the lower limit L by solving numerically Eq. (177|) . for certain values 
of f3 3 and j3§. Note, that a detailed investigation of Eq. f!77|) is rather involved 
and unnecessary for this study However, we have performed computations 
for specific values of f3 3 and /?5, and we present our results in Fig. [TJ In this 
figure he have plotted the function 

H(y) = C 2 (y)-p 3 C 3 (y)-kC s (y) 

versus y, for fixed (3 3 = 5, and several values of (3 5 = —1, 0, 1, 1.6, 1.8, 2, 2.2, 2.5. 
We observe that for /3 5 < 1.6 Eq. (1771) has no real solutions, for 1.6 < < 1.8 
Eq. (1771) has two real solutions, and for 1.8 < (3 5 Eq. (!77j) has one real 
solution. This situation seems to be quite general, as we have performed 
computations for other values of (3 3 , which are not presented in this paper, 
and we observed the same behavior for the function H(y). According to the 
above mentioned results, there is a region of the parameter space for which 
there is no lower limit L, which vanish the Hamiltonian constraint of Eq. 
(177|) . Also, there is a region in which there is an arbitrariness is the choice 
of L, as Eq. (1771) has two distinct solutions for y = L/M. Finally, there is a 
region for which L is determined uniquely by the couplings (3 3 , and fi§ of the 
spherical symmetric object, as Eq. (1771) has one real solution. 

It would be nice if we could determine the range of r by physical and 
geometrical considerations. However, the spherically symmetric solutions 
with zero shift function should satisfy the Hamiltonian constraint. Note that 
the potential in Eq. (1711) diverges for r = 0, so that it is unavoidable to 
introduce a lower limit L in the corresponding integral of Eq. (1751) . As we 
show explicitly above, this lower limit is determined by the free parameters 
of the model, if we try to satisfy the Hamiltonian constraint of Eq. (1751) (or 
of Eq. (173]) ). It seems that the range of acceptable values for r can not be 
determined geometrically, for spherically symmetric solutions with zero shift 
function in Covariant Horava Lifshitz Gravity, unlike the case of standard 
black hole solutions in General Relativity. However, we think that this topic 
needs further investigation. 
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Figure 1: H(y) = c^CMy) — $^Cz{y) — feC^y) versus y, for ^3 = 5 and increasing 
values for ^5, from -1 (upper curve), 0, 1, 1.6, 1.8, 2, 2.2, 2.5 (lower curve). 

5 Conclusions 

We studied the most general case of spherically symmetric vacuum solutions 
in the framework of Covariant Horava Lifshitz Gravity. More specifically, we 
found two classes of spherically symmetric solutions: (i) with nonzero radial 
shift function (see Sec. 14. ip and (ii) with zero radial shift function (see 
Sec. I4.2p . We show that solutions of class (ii) becomes physically relevant if 
we adopt the position of HM for the geometrical role of A in the IR, while 
solutions of class (%) are consistent with observations only if A is independent 
of spacetime geometry. 

In the case of solutions with nonzero radial shift function we would like 
to mention the freedom of the non-dynamical field A(r), although the latter 
has to satisfy constraint equations, as we see in Sec. 14.11 In particular, A(r) 
can be a constant, or it must satisfy a set of constraints which are presented 
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in Sec. 14.1.21 This class of solutions permits an alternative interpretation 
for the field A other than that of HM, as the Newtonian potential for large 
distances is recovered by the radial shift function. In the latter case - situ- 
ation (ii) - one needs to interpret geometrically the U(l) transformation as 
a space time symmetry, in order to recover the Newtonian potential in the 
non-relativistic limit, which in the situation (i) is not necessary. 
Note that, in order to recover the standard Schwarzschild geometry, expressed 
in Painleve-Gullstrand coordinates, one needs to take the solution for which 
A is constant. For the other solutions, where A(r) is not constant and sat- 
isfies the constraints in Sec J4. 1.2] we obtain solutions which obey a modified 
power law for the Newtonian potential, with a leading term 0(r) ~ r~ b and 
b > 1. These solutions may possibly be consistent with solar system tests 
if b is chosen appropriately closely to unity (b ~ 1), but this would still 
need to be checked. However, in Sec J4.1.2l we present a solution, which can 
satisfy solar system tests and corresponds to a nonconstant field A (see Eq. 
( 159]) ). and a spacetime geometry which agree with Schwarzschild metric only 
asymptotically. 

Note that, although the leading order of large distance behavior for the radial 
shift function of A(r) is fixed by the requirement of recovery of Newton Law 
in the large distance limit, for small distances there is an ambiguity in the 
exact shape of A(r), as we see from the examples in Sec. 14.1.21 
Finally, for solutions with non-zero radial component of shift functions, higher 
order spatial curvature corrections in the action do not play a role, and the 
coupling constants in Eq. f )25|) do not appear in this class of solutions. 

In the case of zero radial shift function, in contrast with the previous 
case, higher order curvature terms do not vanish and, as a result, the corre- 
sponding higher order couplings appear in the solution for the auxiliary field 
A. In Sec. I4.2.1l we explain why only the IR part of A possesses a geometric 
interpretation in the IR, while the role of the UV part of A, in the framework 
of HM view, remains ambiguous. 

In addition, this situation is characterized by a difficulty related to the Hamil- 
tonian constraint. In order to satisfy the latter, we had to assume the exis- 
tence of a spherical gap in space (with radius L), centered on the spherical 
compact object. The radius L of the gap can be determined by an algebraic 
equation which can be solved only numerically, as it is shown in Sec. 14.2.21 
Note, that there are regions of the free parameter space for which this equa- 
tion has no solutions, or has more than one solution (in particular two). In 
this sense, for the case of two solutions, there is an arbitrariness in the choice 
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of L. In addition, in this class the auxiliary field A is determined unambigu- 
ously by the corresponding constraint (for details see Sees. 131 and H~2"j) . Also, 
note that this class of solutions can pass solar system tests, failing to impose 
restrictions to Horava couplings. 

The main conclusion of this paper, is that beyond the geometrical approach 
of HM, there is an alternative view for the field A which is also consistent 
with observations. In the latter case the field A is considered to be indepen- 
dent of the spacetime geometry, in contrast to HM theory, and the Newton 
law is reproduced by the nonzero shift function of the solutions of class (i). 
Finally, note that the field A is not completely fixed for solutions of class (i), 
so a study on this topic may be a topic for future investigation. 

Acknowledgments We would like to thank K. Farakos and A. Kehagias for 
numerous valuable discussions. This work is supported by the Royal Society, 
UK. 
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Appendix A: Geometrical interpretation of the 
U(l) symmetry 

We consider the standard metric in the ADM form 

ds 2 = -c 2 lN{t) 2 J dt 2 + 2N l d Xi dt + g^dx'dx 3 , (78) 

in which the lapse function N(t) is assumed to be only a function of time. 
We can promote the auxiliary field A(x,t) as a part of spacetime geometry, 
and simultaneously introduce a new spacetime dependent the lapse function, 
by performing the replacement 

N(t) N(t) - 

c z 



The above metric can be written as 

1 



ds 2 = -c 2 ( N 2 - N Ni + 2AN j dt 2 _| 2X ! (lr l tH-\- fn^l.r'd.r'. (70) 
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where we have dropped the higher order term A 2 /c 4 , as we are interested for 
the nonrelativistic limit c — > oo. 

It is easy to see, in the limit where c — > oo, that the metric f J79|) is 
invariant under a spacetime dependent reparametrization of time, according 
to 

t' = t+^fi-, x' = x, (80) 

c z 

if this transformation is accompanied by the U(l) symmetry 
Ni(x,t) = Ni(x,t) + N 2 ViS 

A'(x,t) = A(x,t)+eN + sN-NN i V i £, (81) 

where 

e{x,t) = 

(Note that, under this U(l) transformation, the fields N(t) and gij(x,t) re- 
main unaltered). Indeed, if we consider the general relativistic metric in the 
i! coordinate 

(N H N' A- 2A'N\ 
N 2 *-± J {dt') 2 + 2N H dxidt' + g VJ dx l dx\ (82) 

where all the fields depend on t', the U(l) transformation of Eqs. (I8ip . 
together with 

W(0 = N (t + %^) = N(t) + (83) 

dt'=(l + ^)dt + ^^ + o(\), (84) 



c 2 J c 2 V c 4 



show that the metrics (1521 and (1791) are equivalent, up to higher order terms 
in 1/c. 



Appendix B: £7(1) symmetry by a minimal sub- 
stitution 

In Ref. [42] is proposed a minimal substitution mechanism, which can be 
used to extend the Gauge symmetry of any Dif f(A4, J 7 ) invariant action. 
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This mechanism is bases on the observation that the following quantities are 
invariant under the U(l) transformation, 

5 a {Ni - NViu) = 0, 5 a (A - a) = (85) 

where a is defined as 

N 

a = v - N 3 Vvj + —VjuVu . (86) 

We consider, a Dif f{M.,J r ) invariant action of the form: 

S[N, Ni, 9ij , i) n ] = S H l[N, N h gij ) + S m [N, N, 9ij , ^ n ), (87) 
in which the first term Shl is the standard HL action 

S HL [N, N, 9ij ] = ^J dtd 3 x^N [KijKV - \K 2 - V] , (88) 

and the second term S m represents the interaction between the external fields 
ip n and the gravitational fields N, N, g^. The fields ip n may be, for exam- 
ple, scalar or vector fields. The action S[N, N, g^, ip n ] can be promoted to 
a manifestly U(l) invariant action S[N,Ni,gij,if) n ,A,u], if we perform the 
replacement Ni — > Ni — NViV, and simultaneously add an extra term which 
depends only on (A — a), according to the equation: 



S[N,Ni,g ih ^ n ,A,v\ = S[N,N i -NV i u,g i3 ,^ n }+ J dtd 3 x^Z{^ n , gij ){A-a) 

(89) 

where Z(ip n , g^), with dimension [Z\ = 2, is the most general operator which 
is invariant under Dif f(M, J 7 ) and respects renormalizability requirements. 
This minimal substitution mechanism give naturally the answer of how one 
can couple matter with the auxiliary fields A and v. 

As shown in Ref. [42J, in the absent of the external field ip n , if we set 
A = 1, the manifestly U(l) invariant action S[N, Ni, g^, 0, A, v\ is identical 
with the action of Covariant HL Gravity, as it is given by Eq. (JHJ) above, or 
equivalently we can write: 

S = S HL [N,N % -NViV,g %j \-^ J ' dtd 3 x^g-(A-a)(R-n) 

-A(#-20)j, (90) 



2 
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Note that, if we set i/j. 



0, we obtain 



Z(Q, 9ij ) = ~(R-n). 



Finally, we observe that this derivation of the Covariant HL action does not 
require A = 1, which implies that U(l) symmetry can not force the value of 
A to be A = 1. However, for A ^ 1, the action of Eq. fl90|) above is modified 
by a term of the form: 
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